Response of a three-species cyclic ecosystem to a short-lived elevation of death rate.

A balanced ecosystem with coexisting constituent species is often perturbed by different natural events that persist only for a finite duration of time. What becomes important is whether, in the aftermath, the ecosystem recovers its balance or not. Here we study the fate of an ecosystem by monitoring the dynamics of a particular species that encounters a sudden increase in death rate. For exploration of the fate of the species, we use Monte-Carlo simulation on a three-species cyclic rock-paper-scissor model. The density of the affected (by perturbation) species is found to drop exponentially immediately after the pulse is applied. In spite of showing this exponential decay as a short-time behavior, there exists a region in parameter space where this species surprisingly remains as a single survivor, wiping out the other two which had not been directly affected by the perturbation. Numerical simulations using stochastic differential equations of the species give consistency to our results.

Chimeras in globally coupled oscillators: A review.

The surprising phenomenon of chimera in an ensemble of identical oscillators is no more strange behavior of network dynamics and reality. By this time, this symmetry breaking self-organized collective dynamics has been established in many networks, a ring of non-locally coupled oscillators, globally coupled networks, a three-dimensional network, and multi-layer networks. A variety of coupling and dynamical models in addition to the phase oscillators has been used for a successful observation of chimera patterns. Experimental verification has also been done using metronomes, pendula, chemical, and opto-electronic systems. The phenomenon has also been shown to appear in small networks, and hence, it is not size-dependent. We present here a brief review of the origin of chimera patterns restricting our discussions to networks of globally coupled identical oscillators only. The history of chimeras in globally coupled oscillators is older than what has been reported in nonlocally coupled phase oscillators much later. We elaborate the story of the origin of chimeras in globally coupled oscillators in a chronological order, within our limitations, and with brief descriptions of the significant contributions, including our personal experiences. We first introduce chimeras in non-locally coupled and other network configurations, in general, and then discuss about globally coupled networks in more detail.

Higher-order interactions in Kuramoto oscillators with inertia.

How do higher-order interactions influence the dynamical landscape of a network of the second-order phase oscillators? We address this question using three coupled Kuramoto phase oscillators with inertia under pairwise and higher-order interactions, in search of various collective states, and new states, if any, that show marginal presence or absence under pairwise interactions. We explore this small network for varying phase lag in the coupling and over a range of negative to positive coupling strength of pairwise as well as higher-order or group interactions. In the extended coupling parameter plane of the network we record several well-known states such as synchronization, frequency chimera states, and rotating waves that appear with distinct boundaries. In the parameter space, we also find states generated by the influence of higher-order interactions: The 2+1 antipodal point and the 2+1 phase-locked states. Our results demonstrate the importantance of the choices of the phase lag and the sign of the higher-order coupling strength for the emergent dynamics of the network. We provide analytical support to our numerical results.

Multimodal distribution of transient time of predator extinction in a three-species food chain.

The transient dynamics capture the time history in the behavior of a system before reaching an attractor. This paper deals with the statistics of transient dynamics in a classic tri-trophic food chain with bistability. The species of the food chain model either coexist or undergo a partial extinction with predator death after a transient time depending upon the initial population density. The distribution of transient time to predator extinction shows interesting patterns of inhomogeneity and anisotropy in the basin of the predator-free state. More precisely, the distribution shows a multimodal character when the initial points are located near a basin boundary and a unimodal character when chosen from a location far away from the boundary. The distribution is also anisotropic because the number of modes depends on the direction of the local of initial points. We define two new metrics, viz., homogeneity index and local isotropic index, to characterize the distinctive features of the distribution. We explain the origin of such multimodal distributions and try to present their ecological implications.

Transition to hyperchaos and rare large-intensity pulses in Zeeman laser.

A discontinuous transition to hyperchaos is observed at discrete critical parameters in the Zeeman laser model for three well known nonlinear sources of instabilities, namely, quasiperiodic breakdown to chaos followed by interior crisis, quasiperiodic intermittency, and Pomeau-Manneville intermittency. Hyperchaos appears with a sudden expansion of the attractor of the system at a critical parameter for each case and it coincides with triggering of occasional and recurrent large-intensity pulses. The transition to hyperchaos from a periodic orbit via Pomeau-Manneville intermittency shows hysteresis at the critical point, while no hysteresis is recorded during the other two processes. The recurrent large-intensity pulses show characteristic features of extremes with their height larger than a threshold and the probability of a rare occurrence. The phenomenon is robust to weak noise although the critical parameter of transition to hyperchaos shifts with noise strength. This phenomenon appears as common in many low dimensional systems as reported earlier by Chowdhury et al. [Phys. Rep. 966, 1-52 (2022)], there the emergent large-intensity events or extreme events dynamics have been recognized simply as chaotic in nature although the temporal dynamics shows occasional large deviations from the original chaotic state in many examples. We need a new metric, in the future, that would be able to classify such significantly different dynamics and distinguish from chaos.

Transition to hyperchaos: Sudden expansion of attractor and intermittent large-amplitude events in dynamical systems.

Hyperchaos is distinguished from chaos by the presence of at least two positive Lyapunov exponents instead of just one in dynamical systems. A general scenario is presented here that shows emergence of hyperchaos with a sudden large expansion of the attractor of continuous dynamical systems at a critical parameter when the temporal dynamics shows intermittent large-amplitude spiking or bursting events. The distribution of local maxima of the temporal dynamics is non-Gaussian with a tail, confirming a rare occurrence of the large-amplitude events. We exemplify our results on the sudden emergence of hyperchaos in three paradigmatic models, namely, a coupled Hindmarsh-Rose model, three coupled Duffing oscillators, and a hyperchaotic model.

Extreme events in a complex network: Interplay between degree distribution and repulsive interaction.

The role of topological heterogeneity in the origin of extreme events in a network is investigated here. The dynamics of the oscillators associated with the nodes are assumed to be identical and influenced by mean-field repulsive interactions. An interplay of topological heterogeneity and the repulsive interaction between the dynamical units of the network triggers extreme events in the nodes when each node succumbs to such events for discretely different ranges of repulsive coupling. A high degree node is vulnerable to weaker repulsive interactions, while a low degree node is susceptible to stronger interactions. As a result, the formation of extreme events changes position with increasing strength of repulsive interaction from high to low degree nodes. Extreme events at any node are identified with the appearance of occasional large-amplitude events (amplitude of the temporal dynamics) that are larger than a threshold height and rare in occurrence, which we confirm by estimating the probability distribution of all events. Extreme events appear at any oscillator near the boundary of transition from rotation to libration at a critical value of the repulsive coupling strength. To explore the phenomenon, a paradigmatic second-order phase model is used to represent the dynamics of the oscillator associated with each node. We make an annealed network approximation to reduce our original model and, thereby, confirm the dual role of the repulsive interaction and the degree of a node in the origin of extreme events in any oscillator associated with a node.

Instabilities in quasiperiodic motion lead to intermittent large-intensity events in Zeeman laser.

We report intermittent large-intensity pulses that originate in Zeeman laser due to instabilities in quasiperiodic motion, one route follows torus-doubling to chaos and another goes via quasiperiodic intermittency in response to variation in system parameters. The quasiperiodic breakdown route to chaos via torus-doubling is well known; however, the laser model shows intermittent large-intensity pulses for parameter variation beyond the chaotic regime. During quasiperiodic intermittency, the temporal evolution of the laser shows intermittent chaotic bursting episodes intermediate to the quasiperiodic motion instead of periodic motion as usually seen during the Pomeau-Manneville intermittency. The intermittent bursting appears as occasional large-intensity events. In particular, this quasiperiodic intermittency has not been given much attention so far from the dynamical system perspective, in general. In both cases, the infrequent and recurrent large events show non-Gaussian probability distribution of event height extended beyond a significant threshold with a decaying probability confirming rare occurrence of large-intensity pulses.

Bistability in a tri-trophic food chain model: Basin stability perspective.

The most important issue of concern in a food chain is the stability of species and their nature of persistence against system parameter changes. For understanding the stable dynamics and their response against parameter perturbation, the local stability analysis is an insufficient tool. A global stability analysis by the conventional techniques seems to supplement some of the shortcomings, however, it becomes more challenging for multistable ecosystems. Either of the techniques fails to provide a complete description of the complexity in dynamics that may evolve in the system, especially, when there is any transition between the stable states. A tri-trophic resource-consumer-predator food chain model has been revisited here that shows bistability and transition to monostability via a border collision that leads to a state of predator extinction. Although earlier studies have partially revealed the dynamics of such transitions, we would like to present additional and precise information by analyzing the system from the perspective of basin stability. By drawing different bifurcation diagrams against three important parameters, using different initial conditions, we identify the range of parameter values within which the stability of the states persists and changes to various complex dynamics. We emphasize the changes in the geometry of the basins of attraction and get a quantitative estimate of the nature of relative changes in the area of the basins (basin stability) during the transitions. Furthermore, we demonstrate the presence of a down-up control, in addition to the conventional bottom-up and top-down control phenomena in the food chain. The application of basin stability in food networks will go a long way for accurate analysis of their dynamics.

Reservoir computing on epidemic spreading: A case study on COVID-19 cases.

A reservoir computing based echo state network (ESN) is used here for the purpose of predicting the spread of a disease. The current infection trends of a disease in some targeted locations are efficiently captured by the ESN when it is fed with the infection data for other locations. The performance of the ESN is first tested with synthetic data generated by numerical simulations of independent uncoupled patches, each governed by the classical susceptible-infected-recovery model for a choice of distributed infection parameters. From a large pool of synthetic data, the ESN predicts the current trend of infection in 5% patches by exploiting the uncorrelated infection trend of 95% patches. The prediction remains consistent for most of the patches for approximately 4 to 5 weeks. The machine's performance is further tested with real data on the current COVID-19 pandemic collected for different countries. We show that our proposed scheme is able to predict the trend of the disease for up to 3 weeks for some targeted locations. An important point is that no detailed information on the epidemiological rate parameters is needed; the success of the machine rather depends on the history of the disease progress represented by the time-evolving data sets of a large number of locations. Finally, we apply a modified version of our proposed scheme for the purpose of future forecasting.

Mitigating long transient time in deterministic systems by resetting.

How long does a trajectory take to reach a stable equilibrium point in the basin of attraction of a dynamical system? This is a question of quite general interest and has stimulated a lot of activities in dynamical and stochastic systems where the metric of this estimation is often known as the transient or first passage time. In nonlinear systems, one often experiences long transients due to their underlying dynamics. We apply resetting or restart, an emerging concept in statistical physics and stochastic process, to mitigate the detrimental effects of prolonged transients in deterministic dynamical systems. We show that resetting the intrinsic dynamics intermittently to a spatial control line that passes through the equilibrium point can dramatically expedite its completion, resulting in a huge reduction in mean transient time and fluctuations around it. Moreover, our study reveals the emergence of an optimal restart time that globally minimizes the mean transient time. We corroborate the results with detailed numerical studies on two canonical setups in deterministic dynamical systems, namely, the Stuart-Landau oscillator and the Lorenz system. The key features-expedition of transient time-are found to be very generic under different resetting strategies. Our analysis opens up a door to control the mean and fluctuations in transient time by unifying the original dynamics with an external stochastic or periodic timer and poses open questions on the optimal way to harness transients in dynamical systems.

Infection spreading and recovery in a square lattice.

We investigate spreading and recovery of disease in a square lattice, and, in particular, emphasize the role of the initial distribution of infected patches in the network on the progression of an endemic and initiation of a recovery process, if any, due to migration of both the susceptible and infected hosts. The disease starts in the lattice with three possible initial distribution patterns of infected and infection-free sites, viz., infected core patches (ICP), infected peripheral patches (IPP), and randomly distributed infected patches (RDIP). Our results show that infection spreads monotonically in the lattice with increasing migration without showing any sign of recovery in the ICP case. In the IPP case, it follows a similar monotonic progression with increasing migration; however, a self-organized healing process starts for higher migration, leading the lattice to full recovery at a critical rate of migration. Encouragingly, for the initial RDIP arrangement, chances of recovery are much higher with a lower rate of critical migration. An eigenvalue-based semianalytical study is made to determine the critical migration rate for realizing a stable infection-free lattice. The initial fraction of infected patches and the force of infection play significant roles in the self-organized recovery. They follow an exponential law, for the RDIP case, that governs the recovery process. For the frustrating case of ICP arrangement, we propose a random rewiring of links in the lattice allowing long-distance migratory paths that effectively initiate a recovery process. Global prevalence of infection thereby declines and progressively improves with the rewiring probability that follows a power law with the critical migration and leads to the birth of emergent infection-free networks.

Routes to extreme events in dynamical systems: Dynamical and statistical characteristics.

Intermittent large amplitude events are seen in the temporal evolution of a state variable of many dynamical systems. Such intermittent large events suddenly start appearing in dynamical systems at a critical value of a system parameter and continues for a range of parameter values. Three important processes of instabilities, namely, interior crisis, Pomeau-Manneville intermittency, and the breakdown of quasiperiodic motion, are most common as observed in many systems that lead to such occasional and rare transitions to large amplitude spiking events. We characterize these occasional large events as extreme events if they are larger than a statistically defined significant height. We present two exemplary systems, a single system and a coupled system, to illustrate how the instabilities work to originate extreme events and they manifest as non-trivial dynamical events. We illustrate the dynamical and statistical properties of such events.

Understanding the origin of extreme events in El Niño southern oscillation.

We investigate a low-dimensional slow-fast model to understand the dynamical origin of El Niño southern oscillation. A close inspection of the system dynamics using several bifurcation plots reveals that a sudden large expansion of the attractor occurs at a critical system parameter via a type of interior crisis. This interior crisis evolves through merging of a cascade of period-doubling and period-adding bifurcations that leads to the origin of occasional amplitude-modulated extremely large events. More categorically, a situation similar to homoclinic chaos arises near the critical point; however, atypical global instability evolves as a channellike structure in phase space of the system that modulates variability of amplitude and return time of the occasional large events and makes a difference from the homoclinic chaos. The slow-fast timescale of the low-dimensional model plays an important role on the onset of occasional extremely large events. Such extreme events are characterized by their heights when they exceed a threshold level measured by a mean-excess function. The probability density of events' height displays multimodal distribution with an upper-bounded tail. We identify the dependence structure of interevent intervals to understand the predictability of return time of such extreme events using autoregressive integrated moving average model and box-plot analysis.

Emergence of Mixed Mode Oscillations in Random Networks of Diverse Excitable Neurons: The Role of Neighbors and Electrical Coupling.

In this paper, we focus on the emergence of diverse neuronal oscillations arising in a mixed population of neurons with different excitability properties. These properties produce mixed mode oscillations (MMOs) characterized by the combination of large amplitudes and alternate subthreshold or small amplitude oscillations. Considering the biophysically plausible, Izhikevich neuron model, we demonstrate that various MMOs, including MMBOs (mixed mode bursting oscillations) and synchronized tonic spiking appear in a randomly connected network of neurons, where a fraction of them is in a quiescent (silent) state and the rest in self-oscillatory (firing) states. We show that MMOs and other patterns of neural activity depend on the number of oscillatory neighbors of quiescent nodes and on electrical coupling strengths. Our results are verified by constructing a reduced-order network model and supported by systematic bifurcation diagrams as well as for a small-world network. Our results suggest that, for weak couplings, MMOs appear due to the de-synchronization of a large number of quiescent neurons in the networks. The quiescent neurons together with the firing neurons produce high frequency oscillations and bursting activity. The overarching goal is to uncover a favorable network architecture and suitable parameter spaces where Izhikevich model neurons generate diverse responses ranging from MMOs to tonic spiking.

Extreme events in a network of heterogeneous Josephson junctions.

We report intermittent large spiking events in a heterogeneous network of forced Josephson junctions under the influence of repulsive interaction. The response of the individual junctions has been inspected instead of the collective response of the ensemble, which reveals the large spiking events in a subpopulation with characteristic features of extreme events (EE). The network splits into three clusters of junctions, one in coherent libration, one in incoherent rotational motion, and another subpopulation originating EE, which resembles a chimeralike pattern. EE migrates spatially from one to another subpopulation of junctions with the repulsive strength. The origin of EE in a subpopulation and chimera pattern is a generic effect of distributed damping parameter and repulsive interaction, which we verify with another network of the Liénard system. EE originates in the subpopulation via a local riddling of in-phase synchronization. The probability density function of event heights confirms the rare occurrence of large events and the return time of EE as expressed by interevent intervals in the subgroup follows a Poisson distribution. The mechanism of the origin of such a unique clustering is explained qualitatively.

Intermittent large deviation of chaotic trajectory in Ikeda map: Signature of extreme events.

We notice signatures of extreme eventslike behavior in a laser based Ikeda map. The trajectory of the system occasionally travels a large distance away from the bounded chaotic region, which appears as intermittent spiking events in the temporal dynamics. The large spiking events satisfy the conditions of extreme events as usually observed in dynamical systems. The probability density function of the large spiking events shows a long-tail distribution consistent with the characteristics of rare events. The interevent intervals obey a Poissonlike distribution. We locate the parameter regions of extreme events in phase diagrams. Furthermore, we study two Ikeda maps to explore how and when extreme events terminate via mutual interaction. A pure diffusion of information exchange is unable to terminate extreme events where synchronous occurrence of extreme events is only possible even for large interaction. On the other hand, a threshold-activated coupling can terminate extreme events above a critical value of mutual interaction.

Asymmetry in initial cluster size favors symmetry in a network of oscillators.

Counterintuitive to the common notion of symmetry breaking, asymmetry favors synchrony in a network of oscillators. Our observations on an ensemble of identical Stuart-Landau systems under a symmetry breaking coupling support our conjecture. As usual, for a complete deterministic and the symmetric choice of initial clusters, a variety of asymptotic states, namely, multicluster oscillation death (1-OD, 3-OD, and m -OD), chimera states, and traveling waves emerge. Alternatively, multiple chimera death (1-CD, 3-CD, and m -CD) and completely synchronous states emerge in the network whenever some randomness is added to the symmetric initial states. However, in both the cases, an increasing asymmetry in the initial cluster size restores symmetry in the network, leading to the most favorable complete synchronization state for a broad range of coupling parameters. We are able to reduce the network model using the mean-field approximation that reproduces the dynamical features of the original network.

Dragon-king-like extreme events in coupled bursting neurons.

We present evidence of extreme events in two Hindmarsh-Rose (HR) bursting neurons mutually interacting via two different coupling configurations: chemical synaptic- and gap junctional-type diffusive coupling. A dragon-king-like probability distribution of the extreme events is seen for combinations of synaptic coupling where small- to medium-size events obey a power law and the larger events that cross an extreme limit are outliers. The extreme events originate due to instability in antiphase synchronization of the coupled systems via two different routes, intermittency and quasiperiodicity routes to complex dynamics for purely excitatory and inhibitory chemical synaptic coupling, respectively. For a mixed type of inhibitory and excitatory chemical synaptic interactions, the intermittency route to extreme events is only seen. Extreme events with our suggested distribution is also seen for gap junctional-type diffusive, but repulsive, coupling where the intermittency route to complexity is found. A simple electronic experiment using two diffusively coupled analog circuits of the HR neuron model, but interacting in a repulsive way, confirms occurrence of the dragon-king-like extreme events.